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Abstract 

It is suggested that a set of positive- and negative-energy oscillations can be reso- 
nantly excited in the inner region of deformed (warped or eccentric) relativistic disks. 
In this paper we examine how a dissipative process affects on this wave excitation 
process. The results show that when the resonant condition in frequency is roughly 
satisfied and thus the oscillations are excited, introduction of a dissipative process 
works so as to decrease the growth rate of the oscillations. When the frequency dif- 
ference of the two oscillations deviates more than a certain amount from that required 
by resonant condition, however, the oscillations are excited by introduction of dissi- 
pative process. This excitation by dissipative process can be understood as a special 
example of the double- diffusive instability. 

Key words: accretion, accretion disks — high-frequency quasi-periodic oscilla- 
tions — oscillations — relativity — resonance — stability — warp — X-rays; stars 



1. Introduction 

Investigation of the origin of high-frequency quasi-periodic oscillations (HF QPOs) ob- 
served in neutron-star and black-hole low-mass X-ray binaries (LMXBs) is one of important 
subjects of astrophysical disks, since it will give clues to understand the structure of the inner- 
most region of relativistic disks, and to estimate spins of the central sources. One of possible 
causes of HF QPOs is disk oscillations in relativistic disks [see Kato (2001) or Kato et al. 
(2008) for review of disk oscillations in relativistic disks]. Related to this subject, Kato (2004, 
2008a, b) suggested that a pair of positive- and negative- energy oscillations can be excited in 
deformed (warped or eccentric) disks by their resonant interaction through disk deformation. 
The presence of this excitation mechanism seems to be confirmed by numerical calculations 
(Ferreira and Ogilvie 2008; Oktariani et al. 2010) and by analytical examinations (Kato et al. 
2011) made by a different method from Kato's original one. 



1 



There remains, however, an important disagreement between analytical results by Kato' 
group (Kato 2004, 2008a, b; Kato et al. 2011) and numerical ones by Ferreira and Ogilvie 
(2008). In Ferreira and Ogilvie's ones, a dissipative process acting on one of oscillations is 
necessary to excite the pair of the oscillations, while in Kato's ones any energy dissipative 
process is unnecessary to excite the set of positive- and negative-energy oscillations. 

Let us consider two oscillation modes the set of whose frequencies and azimuthal wave 
numbers are, respectively, (cji, mi) and (cjz, m 2 ). Disk deformation is denoted by (ud, mo), 
and wd is fixed in the present problem. In the case where the disk deformation is a stationary 
warp, for example, wd = and tuq = 1. Necessary conditions of resonant interaction between 
the above two oscillations through the disk deformation are 

m 2 — m i i mo and u>2 ~ 0J\ ± Wd • (1) 

In Kato's original studies (2004, 2008a, b), he showed that a set of oscillations satisfying condi- 
tions (1) are resonantly excited when their wave energies have opposite signs, assuming that a 
Lindblad resonance occurs in the propagation region of the oscillations. In a subsequent study 
by Kato et al. (2011) (hereafter, Paper I), it is shown that the presence of Lindblad resonance 
in wave propagation region is not essential for excitation. What is essential is opposite signs 
of wave energy of resonantly interacting oscillations with a sufficient accuracy of U2 ~ oo\ ±ud- 
Numerical calculations made by Ferreira and Ogilvie (2008), however, suggests that in addition 
to the above necessary conditions, a dissipative process in one of oscillations is necessary to 
excite the set of oscillations. As the dissipative process, they consider the corotation damping. 
Except for p-mode oscillations with no node in the vertical direction, non-axisymmetric oscil- 
lations are known to be damped by the presence of corotation resonance, even if the corotation 
point (radius) appears in the evanescent region of the oscillations. To mimic this corotation 
damping, they take into account an artificial resistive force acting on oscillations in their nu- 
merical calculations. 

In this paper we take into account a resistive force, as Ferreira and Ogilvie (2008) did, 
into the formulation, and examine how in our formulation the excitation of the set of possitive- 
and negative-energy oscillations is affected by the presence of the dissipative process. The 
results of analyses show that the difference between Kato et al. (Kato 2004, 2008a, b: Kato 
et al. 2011) and Ferreira and Ogilvie (2008) comes from seeing different aspects of the same 
resonant excitation process. 

2. Basic Hydrodynamical Equations and Resistive Damping 

We follow the same procedures as those in Paper I, except that in the present study a 
resistive force acting on a fluid element is additionally considered. This resistive force is assumed 
to work in perturbed states in the direction opposite to the Lagrangian velocity perturbation 
Su. In this case, the hydrodynamical equation describing perturbations over the equilibrium 
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state is 

rn?i li 

-P5u, (2) 



Du m 1^ 

+ + - Vp 



where <5 denotes the Lagrangian variation and D/Dt is the Lagrangian time derivative along 
the perturbed flow u, and (3 is the coefficient of the resistive force and taken to be a positive 
constant. Other notations in equation (2) have their usual meanings. We are interested in 
relativistic disks, but, for simplicity, the Newtonian formulation is adopted. The effects of the 
relativity is mimiced by using the pseudo-Newtonian potential introduced by Paczyhski and 
Wiita (1980). 

Now, we introduce the Lagrangian displacement vector, £, and express equation (2) in 
terms of £. Since 8u is related to the Lagrangian time derivative along the unperturbed flow 
Uo(r), i.e., D^/Dt, by 5u = Do^/Dt, and 5 and D/Dt are commutative (Lynden-Bell and 
Ostriker 1967), we have 

(Du\ D Dl 
5 {-Dt) = Di 6u= Df^ (3) 
Hence, equation (2) is written as 

ii ? + i (v, + iv P )^- fl M. (4) 

A general expression for the right-hand side of equation (4) in terms of the displacement 
vector £ has been obtained by Lynden-Bell and Ostriker (1967) in the case where the perturba- 
tions are adiabatic and have small amplitudes. The results show that po( r ) times the left-hand 
side of equation (4) can be expressed as 

p ||f + 2poK-V)^ + £(£), (5) 

where d/dt is the Eulerian time derivative, po(r) is the density in the unperturbed state, and 
L(£) is a linear Hermitian operator (see in detail Lynden-Bell and Ostriker 1967). Detailed 
expression for L(£) is unnecessary here except that it is a Hermitian. 

Now, we consider the case where perturbations are weakly non-linear. In this case 
equation (4) is written as 



p ^f + 2poK • V)^ + L(t) = -Ppc 



t + PoC{t,t), (6) 



where C(£,£) is the weakly non-linear term. No detailed expression for C(£,£) is given here 
[see Kato (2004, 2008a) for detailed expressions], but an important characteristics of C(£,£) 
is that we have commutative relations (Kato 2008a) for an arbitrary set of rji(r), rj 2 {r), and 
Vs(r), e.g., 

J PoViC(rj 2 ,ri 3 )dV = J p rj 1 C(r l3 ,r l2 )dV = J p r] 3 C(r ll ,r] 2 )dV, (7) 

where the integration is made over the whole volume. As shown later, the presence of these 
commutative relations leads to a simple expression of instability criterion. We suppose that the 
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presence of these commutative relations is a general property of conservative systems beyond 
the assumption of weak nonlinearity. 

Equation (6) is the basic equation to be treated in this paper, which is reduced to the 
one used in Paper 1 in the limit of — 0. Before examining the effects of the /3-term on 
resonant excitation of oscillations, we briefly check whether in the case of C = 0, the /3-term of 
equation (6) really leads to damping of oscillations. To do so, let us take the time-dependence 
of oscillations to be exp(iut). Then, in the case of no damping (i.e., C = 0), equation (6) 
becomes 

-co 2 p $, + 2itup (u • V)£ + L(£) + /3p [iu + (uq • V)]£ = 0. (8) 

This equation is multiplied by the complex conjugate of £ and integrated over the whole volume 
of disks, assuming that po vanishes on disk surface. Taking the imaginary part of the resulting 
equation (using the facts that L and ipo(u ■ V) are Hermitian operators), we have 
1 



-P- (9) 



2 

This shows that the /3-term acts so as to dampen the oscillations as expected 



3. Coupling of Two Oscillations Through Disk Deformation 

As in Paper I, we consider the case where two normal modes of disk oscillations, £i(v,t) 
and $, 2 ( r ^)y resonantly couple through disk deformation, £ D (T",t). Their sets of frequency and 
azimuthal wavenumber are denoted by mi), (0J2, ^2), and (wd, t»d) as mentioned before. 
For example, (ui, mi) means that £i(r,t) is proportional to exp[i(uit — mitp)], where ip is the 
azimuthal coordinate of the cylindrical cordinates (r, </?, z) whose origin is at the disk center 
and the z-axis is in the direction of disk rotation. The resonant conditions are then given by 
equation (1). In general, through the nonlinear coupling term poC(£,£) [see eq. (6)], other 
modes than £ x and £ 2 appear, and their amplitudes as well as those of £ x and £ 2 become 
time-dependent. The purpose here is to examine under what conditions the amplitudes of £ l 
and £ 2 increase with time when they interact each other through a given disk deformation. In 
order to examine this problem, we assume that normal modes of oscillations form a complete 
set, and expand the whole oscillations, £(r,i), resulting from the coupling in the form: 

£(r,t)=A l (t)£ l (r,t)+A 2 (t)Ur,t) + J2Mt)S a (r,t), (10) 

a 

where £ a (v,t), which is proportional to exp[i(u a t — m a ip)], is a normal mode other than ^ l and 
£ 2) and Y^,a represents the summation over such modes. In equation (10) the disk deformation 
is not included, since it is assumed to have a large amplitude so that its time variation can be 
neglected in the present study of excitation of small amplitude oscillations. 
Substitution of equation (10) into equation (6) leads to 



dA 

^ Ipv—^T iui + (uq ■ V) £»+ PiPoAiliui + (uq ■ V)]& 

1,2 



i=l,2 dt 



4 



dA 



i=l,2 z 



(ft 



111) 



In deriving this equation, such small terms as d 2 Ai/dt 2 and fiidAi/dt have been neglected, 
since we are interested in slow secular evolutions of A's due to the /3-term and the coupling 
term. On the right-hand side of equation (11), the coupling terms that are not related to the 
disk deformation are neglected. Furthermore, the coefficients of damping term, /3, of different 
normal modes of oscillations has been distinguished by attaching subscript. This is because, 
as mentioned before, they are introduced to mimic the corotation damping, which depend on 
modes of oscillations. The asterisk * denotes the complex conjugate. 

As discussed in Paper I, the normal modes of oscillations can be regarded approximately 
as a set of orthogonal functions. In particular, in the vertically isothermal disks, orthogonality 
of normal modes approximately holds with the weight function po(r) (Paper I, see also Okazaki 
et al. 1987). That is, for two normal modes, £ 4 and £ with moderately short radial wavelength 
we have 

(M^-> = (M^><%, (12) 

where (po£i£j) is the integration of po£*£j over the whole volume of the disk, the subscripts, i 
and j, being symbols distinguishing different modes. 

Equation (11) is now multipluied by £\ and integrated over the whole volume. Then, 
using the orthogonal relation and the resonant conditions (1), we have an equation describing 
time evolution of Ai(t), which is (for detailed procedures, see Paper I) 

A .E 1 fdA 1 1 

4z— (— - + 

U\ \ dt 2 

= A 2 A D W 12 exp(— iA + t)8 mum2+mD + A 2 Al ) W 12 *exp(— iA-t)6 

TO1,TW2— TOD ' 

(13) 

where 

A + = L0\ — oj 2 — wd and A_ = oj\ — uj 2 + wd- (14) 
The coupling terms W\ 2 and W\ 2 * in equation (13) are time-independent quantities defined by 1 

W 12 = \ ((MI • C(^,e D ))+(MI ■ C(£ D ,£ 2 )))exp(zA + t), (15) 

and 

fa* = \ ((pod ■ c(€ 2 ,e>) \+/pb£ • c(e >) € 2 )\W(tA_t). (i6) 



1 It is noted that W12 and W12* are time-independent quantities, since, for example, (po£* 'C(£2j£d)) an( l 
(po$* ' C(£d>£2)) in equation (15) are proportional to exp(— zA + t). 
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In equation (13), E\ is the wave energy of mode 1 and defined by (e.g., Kato 2001) 



E 1 = ic* 



(17) 



which is a real quantity, since ipo(u • V) is a Hermitian operator. The symbol 5 a ^ in equation 
(13) means that it is the delta-function, i.e., unity when a — b, while zero when a^b. 

Similarly, by multiplying £2 to equation (11) and integrating over the whole volume, we 
have an equation describing time evolution of A 2 (t), which is 

E 2 fdA 2 1 „ „ 

0J2 \ at 2 

= y4iy4 D W r 2 iexp(iA_t)5 m2imi+mD + AiA^W 2 uexp(iA + t)S 

ni2 ,mi —rriY) ; 

(18) 

where W 2 i and W 2 i* are time-independent quantities defiend by 

W 21 = I (U& ■ C(^^ D ))+( Po e 2 ■ C(£ D ,£i)\ W-iA_f), (19) 



W 2U = I ( (M2 • C(£i,&) )+{Pot 2 ■ C(^^i) ) )exp(-zA + t), 



2 

and £2 is the wave energy of mode 2 defined by 



E 2 = -uj 2 



(21) 



U2(poC 2 £ 2 ) ~ «(po^( w o ■ V)£ 2 ) 

In the limit of /?i = /3 2 = 0, equations (13) and (18) are reduced, respectively, to equations (29) 
and (36) in Paper I. 

Here, it is of importance to note that we have the following identical relations: 

W 2l = (W 12 *)\ (22) 

W 2U = (W 12 )*. (23) 

These relations come from the fact that for arbitrary functions, rj 1 , rj 2 , and rj 3 , their order in 
(poVi ' C(r) 2 ,r] 3 )) can be arbitrary changed [see eq. (7) and Kato 2008a]. 

4. Equation Describing Growth Rate 

Equation (13) and (18) show that in the limit of no coupling, A± and A 2 damp, respec- 
tively, by the rate of fii/2 and (3 2 /2, as expected. The next problem is to examine how Ai and 
A 2 evolve with time by the resonant coupling. This is studied in the two cases of m 2 = ni\ +mo 
and m 2 = mi — trd, separately. 

4-1. Case of m 2 = mi + rriD 



New variables defined by 

1 

[ 2 



A 1 (t)=A 1 (t)exp[i-A-t), (24) 



are introduced. Then, from equations (13) and (18), we have 2 , respectively, 

E l fd .1 1 
;— ( — -z-A_ + - 
uj x \dt 2 2 



4i— f4 + 4 A - + ^2)^2 = A^D^2i- (27) 



and 

u 2 \dt 2 2 

Equation (26) and (27) are easily solved by taking A\ and A 2 to be proportional to exp(iat). 
The equation describing a is then found to be 

cr 2 - 4(01 + ft)o- - ^A 2 _ + ilA-CA, - ft) - Iffift - -^-|A D | 2 | W 21 \ 2 = 0, (28) 
2 4 4 4 \brj\hi 2 

where equation (22) has been used. 

In the limit of f3\ = (3 2 = 0, equation (28) is reduced to 

^^^.l', (29) 
and the oscillations are unstable (a 2 < 0) when 

A 2 - + ^rl^ D | 2 |^2i| 2 <0. (30) 

This is the instability condition obtained in Paper I, and (ujiuj 2 / EiE 2 ) < is a necessary 
condition for instability. In another limit of j3i = f3 2 = /3, equation (28) gives 

This equation shows that the term of /3 always works so as to dampen the oscillations, since /3 
appears in the set of a — 2/3/2. 

An interesting and important effect of /3, however, appears when [3 2 ^ fli, which will be 
considered in the next section. 

4-2. Case of m 2 = mi — rrij) 



A 1 (t)=A 1 (t)exp(i-A+t), (32) 



In this case we introduce variables A\ and A 2 defined by 

A 2 (t) = A 2 (t)exJ-i-A + t\ (33) 
Then, from equations (13) and (18), we have, respectively, 

4 *t(s-4 A + + ^'=-™- (34) 

and 

2 In Paper I, we have introduced A\(t) defined by A\(t) = Ai(t)exp(iA-t) , different from equation (24). The 
definition of A\(t) by equation (24) is, however, better, since the resulting equations (26) and (27) are 
symmetric. 
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4i t(l + 4 A + + ^) i2= ^^- (35) 

As in the case of m 2 = mi + mD, the time dependences of A\(t) and A 2 (t) are taken to be 
exp(iat). Then, from equation (34) and (35) we have an equation describing a as 

° 2 - + faa - \a\ + i\A + (/3 2 - ft) - iftft - ^2_\A^\ 2 \W 12 \ 2 = 0, (36) 

where equation (23) has been used. This equation is the same as equation (28) except that A_ 
and W 2 i in equation (28) are now changed to A + and Wi 2 , respectively. 

5. Numerical Calculations of Growth Rate 

Equations (28) and (36) have similar forms. Considering this, we numerically solve the 
following equation: 

a 2 - il/3(l + 8)a - i A 2 + i- A A/3(l -5)- ±(3 2 5 + W = 0, (37) 

where W is an abbreviation of — (cjicu 2 / EiE 2 )\Av\ 2 \W 2 i\ 2 or — (uiu 2 / EiE 2 )\A D \ 2 \Wi 2 \ 2 and 
taken to be positive, since we are interested in the case of (ojiU 2 / EiE 2 ) < 0. Furthermore, /3 
and 5 represent (3 2 and (3i/(3 2 , respectively. In the case of equation (28), A in equation (37) is 
A_, while it is A + in the case of equation (36). 

As mentioned in the previous section, an interesting case is that one of (3i and f3 2 is 
negligibly small compared with the other one. Such a case is really expected. For example, let 
us consider the case where an axisymmetric g-mode oscillation (mi = 0) with one node in the 
vertical direction (rii = 1) resonantly interacts with an one-armed g-mode oscillation [m 2 = 1) 
with two nodes in the vertical direction (n 2 = 2) through a warp (m^ = 1, n-o = I)- 3 111 this 
case, the resonant condition of m 2 = mi + m-Q is satisfied. Furthermore, the resonant condition 
of u>i = w 2 + ud can also be satisfied with w D ~ 0. This is because the propagation region of 
the former oscillation is u 2 < k 2 (k being the epicyclic frequency), while that of the latter one 
is (u — Q) 2 < k 2 , and these two propagation regions overlap in relativistic disks, since n(r) has 
the maximum near r = 4r g (r g being the Schwarzshild radius). The former oscillation has no 
corotation damping = 0), since it is axisymmetric, while the corotation damping occurs 
in the latter oscillation (f3 2 ^ 0). Hence, this case corresponds to 5 = 0, and equation (37) is 
reduced to 

a 2 -i-/3a--A 2 + i-f3A + W = 0. (38) 

The opposite case of f3 2 = and Pij^O may occur, if we consider a different set of oscilla- 
tions. Such cases are formally involved in equation (38) by changing the sign of A. Considering 
these situations, we solve equation (38) for various set of parameters. Dimensionless parameters 
we adopt are A/W 1 / 2 and (3/W 1 ^ 2 . If we obtain the dimensionless frequency, a/W 1 / 2 , whose 

3 This is one of the cases numerically examined by Ferreira and Ogilvic (2008). 
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Fig. 1. Contours of constant (dimensionless) growth rate, — Qa/W 1 ^ 2 , on the (3/W 1 ^ 2 - A/W 1 ^ 2 plane. 
The values labelled to curves denote dimensionless growth rate, — Qa/W 1 ^ 2 . The parameter S is taken to 
be zero. 

imaginary part is negative, the set of oscillations considered are amplified (overstable). 

First, we consider the range of parameters to be examined to understand the charac- 
teristics of equation (38). In the case of ft — 0, equation (38) gives a 2 = A 2 /4 — W, and the 
sign of a 2 changes at A/W 1 ^ 2 = ±2. Hence, the range of —4 < A/W 1 ^ 2 < 4 will be enough to 
examine the basic behavior of equation (38), since both cases of a 2 < (growing oscillation) and 
cr 2 > (constant amplitude oscillations) are involved in this range. The contours of constant 
value of — ^a/W 1 / 2 (i.e., a normalized growth rate of oscillations) in the parameter plane of 
the AjW x l 2 - (3/W 1 / 2 are shown in figure 1 by solving equation (38), i.e., 5 = 0.. 

Figure 1 shows that the growth rate of oscillations is the strongest at the point of A = 
and (3 = on the A/W 1 ^ 2 - (3/W 1 ^ 2 diagram, and decreases as departing from the point on 
the diagram. Next, let us see how the growth rate changes as (3 increases from zero, keeping 
A fixed. Roughly speaking, when lAI/Vr 1 / 2 < 2, the growth rate decreases as (3 increases. In 
the case of IAI/W 1 / 2 > 2, however, the situation is changed. That is, in the latter case we have 
cr 2 = and the oscillations are purely periodic (there is no growth nor damping) in the limit of 
— 0, but the oscillations are amplified as /3 increases from zero. In other words, (3 acts so as 
to excite oscillations. 

In order to see more in detail the effects of (3 on growth rate of oscillations, the /in- 
dependence of growth rate is shown in figure 2 for some values of 6, fixing \A\/W 1 / 2 . The 
figure shows that when 5 = and | A|/jy 1//2 > 2, the (3 terms in equation (38) act so as to excite 
oscillations. If (3/W 1 ^ 2 is sufficiently large, however, the growth rate tends to zero. 

As mentioned before, the simultaneous presence of (3i and /?2 acts so as to dampen the 
oscillations. To demonstrate this numerically, the cases of 5 = 0.3 and 5 = 1.0 are added in 
figure 2, by solving equation (37) [not equation (38)]. 
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Fig. 2. Dimensionless growth rate, — ^sa/W 1 / 2 , as functions of parameter (3/W 1 / 2 for some values of 
A/W 1 / 2 and 5. The thick curves are for 8 — (i.e., f3\ = 0), and the four curves are for AjW 1 ! 2 = 0, 1.0, 
2.0, and 3.0 from the top to the bottom. The thin curves in the middle region of the figure are for 8 = 0.3, 
and the values of A/W 1 / 2 are 0, 1.0, 2.0, and 3.0 from the top to the bottom. The three curves in the 
left-lower corner of the figure are for 8 = 1.0, and the values of A/W 1 / 2 arc 0, 1.0, and 2.0 from the top 
to the bottom (the curve for A/W 1 / 2 = 3.0 cannot be distinguished from that of A/W 1 / 2 = 2.0 on the 
diagram) . 

6. Summary of Growth Rate and Interpretation 

First, let us consider the case of no dissipation (/3 = 5 = 0). In this case the two oscillations 
are independent and have their own frequencies oj\ and u 2 , when their frequencies are far from 
the resonant conditions (1). Here, we examine only the case where the resonant conditions 
are u)\ ~ u 2 + u>d and mi = m 2 + tuq with Wd > 0, in order to avoid unnecessary complication 
considering other cases simultaneously. If |A + | decreases, for example, by a change of ud, 
the frequencies of two oscillation modes, oj\ and lo 2 , move toward oj\ — A + /2 and u 2 + A + /2, 
respectively, so that the exact resonant condition of frequencies is satisfied by these shifted 
frequencies, i.e., (coi — A + /2) = (oj 2 + A + /2) +cjd- 4 This is realized finally at |A + | = 2W 1 ^ 2 . 
If |A + | decreases further, the frequencies of oscillations no more change (frequency rocking), 
but the oscillations become overstable (amplified). This is the results of Paper I. The cause 
of this overstability is energy exchange between positive- and negative-energy oscillations by 
resonance. 

In the case of |A + | > 2W 1 ^ 2 with no dissipation, the oscillations are purely time- 
dependent although frequencies are shifted from their proper ones, U\ and u 2 , as mentioned 
above. Now, we consider what happens if dissipative processes are added (i.e., /3's^ 0) to the 
above case of |A + | > 2W 1 ^ 2 . The results show that the oscillations are amplified if a dissipative 
process works mainly on one of the oscillations (see figure 2). This result can be understood as 

4 For example, we have £ x oc cxp(iu)it)Ai oc [i(o>i — A+ /2)t]Ai oc cxp[z(o;i — A+/2 + a)t), where a = (A+ /2)(1 — 
AW/A^) 1 / 2 and tends to zero as A + /W 1 ' 2 approaches ±2. 
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follows. 

Let us consider the behavior of the oscillation of $, 1 (mode 1). As understood from the 
first paragraph of this section, two forces act on mode 1. One is the restoring force (hereafter, 
force A) which acts so as to make mode 1 oscillate with its proper frequency u\. The other one 
is a force (hereafter, force B) which works so that a fluid element of the oscillation departs from 
the equilibrium state, i.e., works in the direction to make the oscillation unstable. This latter 
force (force B) comes from the resonant coupling with the oscillation of mode 2 [see the term 
of the right-hand side of equation (34)]. In the case of |A + | > 2W 1 ^ 2 , the former force (force 
A) is stronger than the latter one (force B), and mode 1 is a constant amplitude oscillation 
when there is no dissipation. If a dissipative force works on mode 2 (not mode 1), the force 
B working on mode 1 becomes weak during an oscillation, since the resonant term depends on 
the amplitude of mode 2. This means that the net restoring force (i.e., force A minus force B) 
at a phase of oscillation increases compared with that at the phase one cycle before due to the 
decrease of force B. In other words, a fluid element associated with the oscillations returns to 
the equilibrium position with a faster velocity compared with that in one cycle before, and the 
amplitude of oscillations increases with time. If this amplification of mode 1 is stronger than 
the weakening of mode 2 by dissipative process, the whole system of two oscillations of modes 1 
and 2 are amplified. This is the reason why dissipative process amplifies the set of oscillations 
of modes 1 and 2. 

If a dissipative process works also on mode 1 (i.e., f3\ ^ 0), the restoring force acting on 
mode 1 (force A) decreases during one cycle of oscillation of mode 1, and thus the net restoring 
force acting on mode 1 decreases, leading to damping of the whole system of oscillations of 
modes 1 and 2. It is noted that in the above arguments the roles of modes 1 and 2 can be 
changed. 

This excitation mechanism of oscillations is quite similar with double-diffusive instability 
known in other fields of astrophysics and oceanography. For example, overstable convection 
due to double diffusion processes is known in magnetized or rotating media. Let us consider, 
for instance, a stellar atmosphere where the temperature stratification is super-adiabatic but 
convection is suppressed by restoring force of magnetic fields. If thermal conduction works on 
fluid elements and is faster than diffusion process of magnetic field by magnetic diffusivity, the 
net restoring force working on the fluid elements (magnetic restoring force minus buoyancy 
force) increases during an oscillation, since the decrease of buoyancy force is faster than that of 
the magnetic restoring force. Then, we can expect oscillations whose amplitude increases with 
time (overstable convection). This mechanism of overstable convection is first recognized by 
Cowling (1958), although it is first found mathematically by Chandrasekhar (see Chandrasekhar 
1961). Semi-convection also occurs in a chemically stratified medium by thermal conduction 
(Kato 1966), since chemical element diffusion is much slower than thermal diffusion. All of them 
are excitation of overstable oscillations by double-diffusive instability. As recent applications 
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of double- diffusive instability, see, e.g., Charbonnel and Zahn (2007) and Zaussing and Spruit 
(2011) for elements mixing in stars and Latter et al (2010) for protostellar disks. 

7. Discussion 

In order to suggest one of possible origins of high-frequency quasi-periodic oscillations 
observed in NS and BH LMXBs, Kato (2004, 2008a, b) suggested by analytical calculations 
that in deformed (warp or eccentric deformation) disks, a set of positive- and negative-energy 
oscillations can be resonantly excited. Ferreira and Ogilvie (2008) and Oktarian et al. (2010) 
made numerical calculations to examine whether the excitation really exists. Subsequently, 
Kato et al. (2011) (Paper I) made an analytical examination by a method which is different 
from Kato's original one (Kato 2004, 2008a, b). 

The above analytical and numerical investigations now seem to conform that if positive- 
and negative-energy oscillations resonantly interact through disk deformation, they are excited. 
However, there remains an important difference between the results by Kato's group and those 
by Ferreira and Ogilvie. That is, in the results by Ferreira and Ogilvie (2008), an additional 
ingredient is necessary for the excitation. This is a dissipative process acting on one of the 
oscillations. On the other hand, such dissipative process is unnecessary to excite the oscillations 
in analytical and semi-analytical calculations made by Kato and Oktariani et al. 

Our results in this paper seems to resolve this disagreement. Our results in this paper 
show that when frequency deviation from the exact resonant condition is small, i.e., | AI/JV 1 / 2 < 
2, the oscillations are excited without help of any dissipation process. When the deviation is 
large, i.e., IAI/W 1 / 2 > 2, however, a dissipative process ((3 ^ 0) is necessary to excite the 
oscillations. In this context, we think that the case examined by Ferreira and Ogilvie (2008) 
corresponds to the case of IAI/W 1 / 2 > 2. The following seems to also support this conjucture. 
In the numerical results by Ferreira and Ogilvie (2008), the growth rate increases from zero 
with increase of the damping factor /3, but saturates at a certain level (see figure 4 of their 
paper). This is qualitatively consistent with our result shown in figure 2. 5 

It is noted that the analyses made by Kato (2004, 2008a, b) and those made by Paper I 
and this paper are rather different. In the former studies, the presence of a Lindblad resonance 
in the propagation region of oscillations is assumed and the condition of excitation of oscillations 
are derived by considering the behavior of oscillations around the resonance. In the present 
study as well as in Paper I, such assumption is not involved in the analyses. In this sense, 
the latter analyses are more general. However, it should be rememberd that the magnitude of 
growth rate depends on the magnitude of the coupling terms such as |Wi2| 2 and |H^2i| 2 ; and 
they are large when a Lindblad resonance appears in the propagation region of the oscillations. 

5 In our results the growth rate decreases and tends to zero, when the damping factor increases further. In 
figure 4 of Ferreira and Ogilvie (2008), such trend seems to be also present. 
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It is also noted that for the present resonant excitation mechanism to work efficiently, 
the propagation regions of the positive- and negative-energy oscillations must overlap in the 
radial direction. Otherwise, the terms of resonant couplings are small and there is practically 
little excitation. In nonself-gravitating disks, the overlapping of the propagation regions is easy 
to occur in relativistic disks. 

In our excitation mechanism of disk oscillations, a disk deformation (e.g., warp or eccen- 
tric deformation) is necessary. What kinds of origins of such disk deformation are conceivable? 
There are some possibilities. For example, the angular momentum axis of disks will not be 
always expected to align with the spin axis of the central source. Then, the disks are tilted 
and have precession. In such disks, excitation of a set of disk oscillations with positive- and 
negative-energies are expected. It is noted that the direct cause of this excitation is not the 
spin of the central source, but a disk deformation. As some other possible causes of disk de- 
formation, the disk instability due to irradiation from central sources (Pringle 1992) will be 
conceivable. The Papaloizou-Pringle instability also makes the inner region of disks deformed 
at the phase where the inner tori are formed (Machida and Matsumoto 2008). 

Finally, let us briefly mention whether our excitation mechanism can describe the spin- 
frequency relation observed in QPOs of millisecond pulsars. When a burst occurs on the surface 
of spinning central neutron stars, a non-axisymmetric pattern rotating with the spin frequency, 
u s , will be induced on the disks, i.e., wd = u s . In such cases, excitation of oscillations which 
satisfy U\ — u 2 ~ uj s is expected, and QPOs with such frequency relation are really observed 
(e.g., 4U 1702-43). In some other millisecond pulsars, however, the observed frequency relation 
is u 2 — U\ ~w s /2 (e.g., SAX J1808. 4-3658). At a glance, it seems to be difficult to describe the 
latter by the present resonant model, although more considerations will be worthwhile. 

The author thanks Atsuo T. Okazaki for discussions at the initial phase of this study. 
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